Detailed electronic structure studies on superconducting MgB2 

and related compounds 



P.Ravindran, P.Vajeeston, R.Vidya, A.Kjekshus and H.Fjellvag 

^ Department of Chemistry, University of Oslo, Box 1033, Blindern, N-0315, Oslo, Norway 

(February 1, 2008) 



Abstract 

Our recent electronic structure studies on series of transition metal diborides 
indicated that the electron phonon coupling constant is much smaller in these 
materials than in superconducting intermetallics. However experimental stud- 
ies recently show an exceptionally large superconducting transition tempera- 
ture of 40 K in MgB2 . In order to understand the unexpected superconducting 
behavior of this compound we have made electronic structure calculations for 
MgB2 and closely related systems. Our calculated Debye temperature from 
the elastic properties indicate that the average phonon frequency is very large 
in MgB2 compared with other superconducting intermetallics and the excep- 
tionally high Tc in this material can be explained through BCS mechanism 
only if phonon softening occurs or the phonon modes are highly anisotropic. 
We identified a doubly-degenerate quasi-two dimensional key-energy band in 
the vicinity of Ep along F-A direction of BZ (having equal amount of B px and 
Py character) which play an important role in deciding the superconducting 
behavior of this material. Based on this result, we have searched for similar 
kinds of electronic feature in a series of isoelectronic compounds such as BeB2 , 
CaB2, SrB2, LiBC and MgB2C2 and found that MgB2C2 is one potential ma- 
terial from the superconductivity point of view. We have also investigated 
closely related compound MgB4 and found that its Ep is lying in a pseudo- 
gap with a negligibly small density of states at Ep which is not favorable 
for superconductivity. There are contradictory experimental results regard- 
ing the anisotropy in the elastic properties of MgB2 ranging from isotropic, 
moderately anisotropic to highly anisotropic. In order to settle this issue we 
have calculated the single crystal elastic constants for MgB2 by the accurate 
full-potential method and derived the directional dependent linear compress- 
ibility. Young's modulus, shear modulus and relevant elastic properties from 
these results. We have observed large anisotropy in the elastic properties 
consistent with recent high-pressure measurements. Our calculated polar- 
ized optical dielectric tensor shows highly anisotropic behavior even though 
it possesses isotropic transport property. MgB2 possesses a mixed bonding 
character and this has been verified from density of states, charge density and 
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crystal orbital Hamiltonian population analyses. 



Typeset using REVT^ 



I. INTRODUCTION 



The recent discovery of superconductivity with high Tj^ in MgB2 has initiated large ac- 
tivity in experimental as well as theoretical studies. NbsGe is long been the record holder 
of the highest among the intermetallic superconductors. The recently found YPd2B2C 
touches the same Tc as NbsGe. The work of Bednorz and Miilleii in 1986 on the basis of 
copper oxides started the discovery of a rapidly increasing number of high-T^ superconduc- 
tors, among which the current record is held by HgBa2Ca2Cu308-5 with Tc~ 164 K under 
pressure.! During the past decade, remarkable progress in basic research and technological 
applications has been made on the high-T^ cuprate superconductors. However the complex 
crystal structures in combination with the multicomponent nature of the materials involved 
hinder the full understanding of microscopic origin of high-Tc superconductivity. Hence it 
may be beneficial to study the properties of a simple compound like MgB2 in detail, which 
does not only have high T^, but also takes a simple crystal structure with sp electrons 
(involved in the superconducting process) that are easy to handle theoretically. 

In cuprates it is now generally believed that the conductivity takes place in the Cu02 
planes which, therefore, are essential to the high-Tc superconductors. The observation of 
superconductivity with high Tc in MgB2 releases the question of whether superconductivity 
with even higher Tc can be found in intermetallics which do not comprise the characteristic 
Cu02 planes. Several properties of MgB2 appear closely related to high Tc superconducting 
cuprates: a low electron density of states, a layered structural character and the presence 
of rather light atoms (like oxygen in cuprates) facilitate high phonon frequencies. Further, 
although the electronic structure of this material has three-dimensional character, the B-B 
a bands derived from B p^^y electrons (believed to be important for the superconductivity) 
reflect two-dimensional character. Like the cuprates, MgB2 appears to exhibit hole con- 
ductivity, as evidenced from theoretical considerationJi^ and experimental Hall coefficient 
measurement .i Another interesting aspect of the MgB2 structure is that it has negatively 
charged honeycomb- shaped B planes which are reminicent of highly negatively charged Cu-0 
planes in the high-Tc cuprates. 

It is often believed that the d electrons play an important role for the superconducting 
behaviour of intermetallic compounds. So the experimental search for new superconductors 
has to a large extent been focussed on transition metal compounds. The back-ground for 
this is that transition metal compounds usually possess large density of states (DOS) at the 
Fermi level (one of the ingredients required for high Tc) than the main-group metal {sp) 
compounds. The BCS superconductors usually possess larger DOS at the Fermi level [Ep) 
and the high Tc in conventional superconductors is related to the large N{Ef) values as 
well as the strong coupling of selected phonon modes to the electronic system. Further, 
relatively high Tc in superconducting intermetallics are believed to be coupled to a van 
Hove-like peak at Ep in the DOS curve, as predicted by various bandstructure calculations 
and indicated experimentally for Ni-site substituted rare-earth nickel borocarbides.i How- 
ever in MgB2 there is no such peak-like feature at the vicinity of Ep in the DOS profile 
and hence one must search for some other origin of the superconductivity. MgB2 is a sp 
metal and hence its value of N{Ef) is small compared with that of superconducting tran- 
sition metal compounds. Further, the presence of light borons in MgB2 similar to oxygens 
in high-Tc cuprates, indicates that some exotic mechanism is involved in the superconduc- 
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tivity. On the other hand, theoretical studied show that the quasi two-dimensional B a 
bands are strongly coupled with the phonon modes which would be consistent with 
the BCS mechanism. The experimentally observed isotope effect^ in MgB2 also indicates 
phonon-mediated superconductivity. But^s discussed by Baskaran,llll the absence of Hebel- 
Slichter peak in the NMR relaxation ,0^0 a temperature dependent peak around 17meV 
in the energvresolved neutron scattering,lli first order metal-to-metal transition on Al or C 
substitution]!!^ apparently anomalous temperature dependences of the Hall coefiicientlii 
and London penetration depthS and the need of small /x* value to explain the experimen- 
tally observed high T^. by the BCS theory0 indicate that the mechanism differs significantly 
from that of the BCS theory. 

The discovery of new classes of quaternary intermetallic superconductors; the 
borocarbideiB and the boronitrideil series, with relatively high Tc (up to 23 K for 
YPd2B2Cii) has encouraged the search for superconductivity in materials possessing light 
atoms such as B, C, N, H etc. Superconductivity with the T^ ranging 2—4 K has been 
observed for ternary transition metal diborides (YRe2B2,LuB2C2, YB2C2)I^@ and it is also 
observed in magnetic rare-earth (RE) rhodium borides (RERh4B4c3). Very recently, the 
layer-structure compounds /3-ZrNCl and /5-HfNCl have been found to be superconducting 
upon Li intercalation with = 12.5 and 25.5 K, respectively.il Another group of materials 
maximum T, = 9.97 K reported^ for Y2C2I2) is RE carbide halide superconductors with 
C— C pairs located in octahedrally coordinated voids of close-packed RE atoms. The recent 
discovery of superconductivity0 at :^14K in the AlB2-type phase of CaSi2 (the highest Tc 
ever obtained among the silicides) indicates that the AlB2-type structure may be favorable 
for superconductivity. 

There are several mechanisms proposed for high superconductivity in MgB2 and the 
like. One is based on band structure findings,@ill which suggests that the superconducting 
state results from strong electron-phonon interaction and high phonon frequency associated 
with the light boron atom. This is supported by the recent observation of relatively large 
boron isotope effect on Tc.0 A recent high-pressure studylil shows that Tc decreases with 
pressure at a rate of —1.11 K/GPa which is consistent with the BCS framework. It has been 
widely speculated that the low mass of boron is conducive for the occurrence of high phonon 
frequencies and consequently for high T^. Another mechanism is called the "universal" by 
Hirschli, and this conjectures that the superconductivity in MgB2 (similar to that in cuprate 
superconductors) is driven by the pairing of the heavily dressed holes in bands that are almost 
full to gain enough kinetic energy to overcome the Coulomb repulsion. A positive pressure 
effect on has also been predicted by Hirsch when the pressure reduces the interatomic 
B— B distance. An and PickettS maintain that the B a bands are playing an important role 
in the superconductivity of MgB2, and that the B in-plane E2g phonon mode is strongly 
coupled to this band. Contradictory to the above viewpoint, BaskaranO concluded from 
the resonance valence bond (RVB) theory that the two dimensional a bands do not play 
a crucial role in establishing high superconductivity and pir band might interfere with 
superconductivity. 

MgB2 possesses a hexagonal crystal structure and hence one can expect anisotropy in 
the physical properties of this material. Despite the strongly anisotropic layered hexagonal 
structure, its electronic properties indicate three dimensionality and the calculated Fermi 
velocity in the plane and perpendicular to the plane are almost the same. The theory for 
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hole superconductivityi suggests that a decrease in the B-B interatomic distance should 
increase Tc. The RVB theoryS suggests that a increase in chemical pressure along the c 
axis should decrease and ultimately convert the material to a normal metallic state. In 
order to test these possiblities it is important to know the anisotropy in mechanical prop- 
erties. From high-pressure total energy studies Loa and SyassenE^ concluded that MgB2 
has isotropic compressibility. From a high-resolution x-ray powder diffraction high-pressure 
study along with the density functional calculation Vogt et a/.il concluded that MgB2 pos- 
sesses nearly isotropic mechanical behavior. From isothermal compressibility measurements 
by synchrotron x-ray diffraction Prasside^ concluded that MgB2 is a stiff tightly-packed 
incompressible solid with only moderate bonding anisotropy between inter- and intralayer 
directions. However, Jorgensen et al!^ found unusually large anisotropy in thermal expan- 
sion and compressibility from the neutron diffraction measurements. So, it is interesting 
to calculate the single- crystal elastic constants of MgB2 to identify the exact nature of the 
anisotropy. 

Several attempts have been made to enhance by substitution of such as Al,lli Be,@, 
and Li,ii for Mg and for B but no practical progress has hitherto been obtained. 
Moreover, the role of Mg and B site substitution on the electronic structure of MgB2 has 
been studied theoretically.ii So, it is interesting to search for compounds with an electronic 
structure similar to that of MgB2. A systematic investigation of this and related materials 
is a way to learn the mechanism of superconductivity in this novel material and to identify 
materials with high Tc. Knowledge of electronic structure, DOS, Debye temperature and re- 
lated properties are important for assessing the mechanism and nature of superconductivity. 
In a search for superconducting diborides, BeB2 is a promising candidate since the lighter 
Be may help to provide larger phonon frequencies and hence increase T^. If the electron 
per atom ratio is important for the superconductivity in this class one should also consider 
CaB2 and SrB2 as interesting candidates. If it is the combination of Mg and B which brings 
the superconductivity in MgB2, one has to consider MgB4 also. If the number of electrons 
in the B layers is a key for superconductivity in MgB2, related layer-structured materials 
such as LiBC, MgB2C2 should be paid attention. For these reasons we have made detailed 
electronic structure studies for the above mentioned compounds. 

The rest of the paper is organized as follows. The structural aspects and the compu- 
tational details about the calculations of the electronic structure, optical spectra and the 
elastic constants are given in in Sec.|l[ In Sec. |T| we have analyzed the bonding behav- 
ior of MgB2 using the orbital and site projected DOS, crystal orbital overlap Hamiltonian 
population (COHP), charge density analysis etc. The electronic band structure of MgB2 is 
calculated and compared with that of closely related systems and also analyzed the possible 
connection between the electronic structure and the superconductivity. The elastic and the 
optical anisotropy of this material are calculated and compared with available experimental 
results. Finally we summarize the important findings of the present study in Sec. [V[ 



II. STRUCTURAL ASPECTS AND COMPUTATIONAL DETAILS 
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A. Crystal structure details 

MgB2 (Fig.|I]a) has AlB2-type structured with space group P6/mmm and lattice param- 
eters, a = 3.084 Aand c = 3.522 A. It is a simple hexagonal lattice of close-packed Mg layers 
alternating with graphite- like B layers, viz. B atoms arranged at the corners of a hexagon 
with three nearest neighbor B atoms in each plane. The Mg atoms are located at the center 
of the B hexagon, midway between adjacent B layers. 

MgB2C2 crystallizes in an orthorhombic structure,^ space group Cmca with a = 10.92, 
b = 9.46 and c = 7.45 A. The structure of MgB2C2 (Fig.|l|b) contains graphite-like but 
slightly puckered boron-carbon layers whose charge is counter-balanced by Mg^+ cations. 
The mutual coordination of boron and carbon consists of five atoms of the other kind, three 
of which being located in the same and two in adjacent layers. Each of Mg is coordinated 
by six B and six C atoms arranged at the corners of a slightly distorted hexagonal prism. 
The B-C distances within the layers range from 1.562 to 1.595 A. 

LiBC (Fig. |1]3) crystallizesil in a hexagonal primitive lattice with space group PQ^/mmc. 
The lattice parameters are a = 2.752 and c = 7.058 A. The B and C atoms form a planar 
so-called heterographite layer. The interlayer regions are filled by Li. The B-C distance of 
1.589 Ain LiBC is comparable with that in MgB2C2. 

In MgB4, the B atoms form interconnected pentagonal pyramids with the Mg atoms 
located in channels running parallel to the c axis. The Mg atoms form zig-zae chains. MgB4 
is orthorhombic, space group Pnma, with a = 5.46, b = 7.47 and c = 4.42 AB^. The average 
B-B distance in the pentagonal pyramid is 1.787 A. Whenever possible we have used the 
experimental lattice paramenters for our calculation, whereas for BeB2, CaB2 and SrB2 we 
have used the optimized structural parameters obtained from total energy minimization. 

B. Computation details for the FPLAPW calculations 

These investigations are based on ah initio electronic structure calculations derived from 
density-functional theory. For the screened plasma frequency and the orbital projected DOS 
calculations we have applied the full-potential linearized-augmented plane wave (FPLAPW) 
methodB in a scalar-relativistic version without spin-orbit (SO) coupling. In the calculation 
we have used atomic sphere radii 1.8 and 1.5 a.u. for Mg and B, respectively. The charge 
density and the potentials are expanded into lattice harmonics up to £ = 6 inside the spheres 
and into a Fourier series in the interstitial region. The initial basis set included 3s, 3p and 3d 
valence and 2s, 2p semicore functions at the Mg site, 2s, 2p and 3d valence functions for the B 
site. The set of basis functions was supplemented with local orbitals for additional flexibility 
in representing the semicore states and for relaxing the linearization errors generally. The 
effects of exchange and correlation are treated within the generalized-gradient-corrected 
local-density approximation using the parameterization scheme of Perdew et al^ To ensure 
convergence for the Brillouin zone (BZ) integration, 320 k-points in the irreducible wedge 
of the first BZ of the hexagonal lattice for MgB2 were used. Self-consistency was achieved 
by demanding the convergence of the total energy to be smaller than 10~^Ry/cell. This 
corresponds to a convergence of the charge below 10"'^ electrons/atom. For BeB2, CaB2 and 
SrB2 we have made the structural optimization with the similar procedure. 
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C. Computation details for the TB-LMTO calculations 

To calculate the electronic ground state properties of MgB4, MgB2C2 and LiBC we used 
the TB-LMTO method of Andersen.lil The von Barth-Hedin parameterization is used for 
the exchange correlation potential within the local density approximation. In the present 
calculation, we used atomic sphere approximation. The calculations are semi-relativistic, 
i.e. except spin-orbit coupling, all other relativistic effects are included, taking also into ac- 
count combined correction terms. BZ k-point integrations are made using the tetrahedron 
method on a grid of 405 (MgB4), 365 (MgB2C2) and 549 (LiBC) k-points in the irreducible 
part of BZ. In order to have more insight into the chemical bonding, we have also evalu- 
ated the crystal orbital Hamiltonian population (COIIP)EE in addition to the regular band 
structure calculations. COHP is the density of states weighted by the corresponding Hamil- 
tonial matrix elements, a positive sign of which indicating bonding character and negative 
antibonding character. 

D. Computational details of FPLMTO calculations 

The full-potential LMTO calculationdli presented in this paper are all electron, and no 
shape approximation to the charge density or potential has been used. The base geometry 
in this computational method consists of muffin-tin and interstitial parts. The basis set is 
comprised of augmented linear muffin-tin orbitals.0 Inside the muffin-tin spheres the basis 
functions, charge density and potential are expanded in symmetry adapted spherical har- 
monic functions together with a radial function. Fourier series is used in the interstitial 
regions. In the present calculations the spherical-harmonic expansion of the charge density, 
potential and basis functions were carried out up to ^ = 6. The tails of the basis functions 
outside their parent spheres are linear combinations of Hankel or Neumann functions de- 
pending on the sign of the kinetic energy of the basis function in the interstitial regions. 
For the core-charge density, the Dirac equation is solved self-consistently, i.e., no frozen core 
approximation is used. The calculations are based on the generalized-gradient-corrected- 
density-functional theory as proposed by Perdew et a/.ii SO term is included directly in 
the Hamiltonian matrix elements for the part inside the muffin-tin spheres. Moreover, the 
present calculations make use of a so-called multi basis, to ensure a well converged wave 
function. This means that we use different Hankel or Neuman functions each attaching to 
its own radial function. This is important to obtain a reliable description of the higher lying 
unoccupied states, especially for the optial property studies. For our elastic properties study 
we have used 192 k points and for the optical property studies 624 k points in the irreducible 
part of BZ. 

E. Calculation of optical properties 

Once the energies ekn and functions |kn) for the n bands are obtained self consistently, 
the interband contribution to the imaginary part of the dielectric functions €2(00) can be 
calculated by summing the transitions from occupied to unoccupied states (with fixed k 
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vector) over BZ, weighted with the appropriate matrix element for the probabihty of the 
transition. To be specific, the components of €2(00) are given by 



where {px,Py,Pz)=P is the momentum operator and /kn is the Fermi distribution. The eval- 
uation of matrix elements in Eq. (|l]) is done over the muffin-tin and interstitial regions 
separately. Further details about the evaluation of matrix elements are given elsewhere." 
For hexagonal structure of MgB2 the dielectric function is a tensor. By an appropriate 
choice of the principal axes we can diagonalized it and restrict our considerations to the 
diagonal matrix elements. We have calculated the two components E\\a and E\\c of the 
dielectric constants corresponding to the electric field parallel to the crystallographic axes a 
and c, respectively. These calculations yield the unbroadened functions. To reproduce the 
experimental conditions correctly, it is necessary to broaden the calculated spectra. The 
exact form of the broadening function is unknown, although comparison with measurements 
suggests that the broadening usually increases with increasing excitation energy. Also the 
instrumental resolution smears out many fine features. To simulate these effects the lifetime 
broadening was simulated by convoluting the absorptive part of the dielectric function with 
a Lorentzian, whose full width at half maximum (FWHM) is equal to 0.01 (^tu)^ eV. The ex- 
perimental resolution was simulated by broadening the final spectra with a Gaussian, where 
FWHM is equal to 0.02 eV. For metals, the intraband contribution to the optical dielectric 
tensor influence to the lower energy part of the spectra. This has been calculated using the 
screened plasma frequency obtained from Fermi surface integration with the description in 



The hexagonal phase of MgB2 has two lattice parameters a and c with Bravais lattice 
vectors in matrix form 



The FPLMTOeHI method allows total-energy calculations to be done for arbitrary crystal 
structures. We can therefore apply small strains to the equilibrium lattice, then determine 
the resulting change in the total energy, and from this information deduce the elastic con- 
stants. The elastic constants are identified as proportional to the second-order coefficient 
in a polynomial fit of the total energy as a function of the distortion parameter 5.11 We 
determine linear combinations of the elastic constants by straining the lattice vectors R 
according to the relation R' = RD. Here R' is a matrix containing the components of the 
distorted lattice vectors and D the symmetric distortion matrix, which contains the strain 
components. We shall consider only small lattice distortions in order to remain within the 
elastic limit of the crystal. In the following we shall briefiy list the relevant formulas used to 




(1) 



Ref. 0. 



F. Calculation of Elastic properties 
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obtain the elastic constants for hexagonal crystals. The internal energy of a crystal under 
strain, 6, can be Taylor expanded in powers of the strain tensor with respect to the initial 
internal energy of the unstrained crystal in the following way, 



E{V, 5) = E{Vo, 0) + T/o(E ^i^i^i + 1/2 E ^iA^i^jQ 



(2) 



The volume of the unstrained system is denoted Vo, E(Vo,0) being the corresponding total 
energy. In the equation above, is an element in the stress tensor. 

Since we have five independent elastic constants, wc need five different strains to deter- 
mine them. The five distortions used in the present investigation are described below. The 
first distortion 

fl+6 \ 
Di = 1 + S 

V 1 + S J 

gives compression or expansion to the system. This preserves the symmetry but changes the 
volume. The strain energy associated with this distortion is 



E(V, 5) = E{Vo, 0) + Vo[{n + T2 + rs)5 + -(2cn + 2ci2 + 4ci3 + caa)^'] 



The second distortion 



Do = 



V 





+ 

(1 + 5)2/3 J 



gives volume and symmetry conserving variation of c/a. The energy associated with this 
distortion is 

E{V, d) = E{Vo, 0) + Voiin + r2 + T^)d + ^(cii + cu - 4ci3 + 2033)^'] 

y 



The strain matrix 



/ 1+5 

' (1-52)1/3 



\ 










1-5 



(1-52)1/3 








(1-52)V^ / 



distorts the basal plane by elongation along a and compression along b in such a way that 
the volume is conserved. The energy associated with this distortion is 

E(V, 5) = E(Vo, 0) + Vo[(n - r2)5 + (cn - cu)5^)] 

The elastic constant C55 can be determined by the distortion of the lattice using the 
volume conserving triclinic distortion 
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Da 



( (1-52)1/3 






1 



(1-52) 





1/3 \ 



(1-52)1/3 

\ (1-52)1/3 ^ (1-52)1/3 / 




The energy change associated with this distortion is 

E{V,5) = E(K,,0) + K)[r55+ (2c55)5'] 

The fifth strain 



involves stretching of the c axis while keeping other axes unchanged. Hence, the hexagonal 
symmetry is preserved but volume is changed. The energy change associated with this strain 
can be written as 

E{y,5) = E{y,,Q) + V,[r^5 + C-f)5^] 

The elastic constant C33 can be directly obtained from the above relation. By solving the 
linear equations given above we have obtained all the five elastic constants. From pressure 
dependent lattice parameter measurements it is easy to obtain the bulk modulus along the 
crystallographic axes. Also to quantify the mechanical anisotropy of MgB2 it is important to 
calculate the bulk modulus along the axes. For hexagonal crystals the bulk modulus along 
a (Ba) and c (Be) are defined as 



Ba 



dP 

da 



A 



1 + a 



and 



dP Ba 
Be = c— = — 
dc a 

where A = 2(cii + C12) + 4ci3a + 0330?^ and 

Cll + C12 — 2Ci3 



a 



C33 - Ci3 



The calculated bulk modulus along the crystallographic axes obtained from these relations 
are compared with the experimental results in Sec.|ITT|. 

III. RESULTS AND DISCUSSION 



A. Electronic structure 

The calculated electronic band structure of MgB2 is given in Fig.|^. The interesting 
feature of this band structure is that a doubly degenerate, nearly fiat, bands are present 
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just above the Ep in the F— A direction in Fig.^ and cut the Ep along the K— F direction. 
These bands give rise to nearly cylindrical, hole-like Fermi surfaces around the F point ,@ 
indicating that the transport properties are dominated by the hole carriers in the plane where 
B atoms exist. These bands are incompletely filled bonding a bands with predominantly 
boron p^^y character. The pz bands (they are mainly in the unoccupied state and having 
finite contribution along M-F direction at the VB in Fig. ^ are derived from the intralayer tt 
bonding orbitals which also have interlayer couplings between adjacent atomic orbitals in the 
c direction. Our earlier studylil on superconducting LaaX (X = Al,Ga, In, Tl) compounds 
show that the presence of a fiat band in the vicinity of Ep gives large (For Laain addition 
of carbon gives a stable LaaInC compound for which Ep is brought to the pseudogap and 
results in a non-superconducting state.) The fiat band feature is also present in the recently 
discovered superconducting YNi2B2C compound. As a working hypothesis, we believe that 
this fiat-band feature plays an important role for the superconductivity in MgB2. It is 
worth to recall that the calculatedc^ electron-phonon interaction strength also shows a large 
value in the F-A direction where the fiat-band feature is seen. Moreover, zone boundary 
phonon calculations show that this band feature is very sensitive to the E2g mode (B-bond 
stretching). We observed that the top of this fiat-band feature is around 0.54 eV above Ep. 
Our calculations show that addition of around 0.32 electron to MgB2 will bring Ep to the 
top of this energy band (assuming rigid band filling). Thus, if electrons are responsible for 
the superconductivity one can expect enhancement of on electron doping. Slusky et a/.lli 
have reported the role of electron doping on the superconductivity in Mgi_2:Ali,.B2 phase 
for which it is found that the Tc drops smoothly up to x = 0.1 and beyond x = 0.25 the 
superconductivity is completely destroyed. This indicates that holes are responsible for the 
superconductivity in MgB2 an observation consistent with Hall effect measurements.! The 
electron per atom ratio for TiB2 is similar to that for MgB2 and hence it seems worthwhile 
to make comparison between these phases. The electronic structure of TiB2@ shows that 
Ep is located in a pseudogap on the DOS curve and hence the estimated electron-phonon 
coupling constant is much smaller than for superconducting materials. This observation is 
consistent with the experimental findings in the sense that no superconductivity is detected 
for TiB2 even below 1 K.il 

The role of Mg on the band structure of MgB2 can be elucidated by completely removing 
the Mg atoms from the lattice and repeating the calculations for a hypothetical structure 
with only B atoms (note: using the lattice parameters for MgB2). The calculated electronic 
structure for B network alone is given in Fig.|^ along with that of MgB2. The striking 
difference between the two cases is in the position of the flat band in the F-A direction. 
This flat band feature for the B network alone is strongly two dimensional (viz. very little 
dispersion along F-A). This is ca. 1.8 eV above Ep owing to the lower number of electrons 
in the B network compared with MgB2. This suggests that the B-B a bonds are primarily 
responsible for the superconductivity in MgB2. Interestingly the bands are not deformed 
appreciably when we remove Mg from MgB2 indicating that the electrons from Mg atoms 
mainly give a shift in Ep almost like rigid-band fllling. This view point is conflrmed also 
from the density of states study (see Fig.,^ which shows that the topology of the DOS 
proflle with and without Mg atoms in MgB2 are almost the same. But a shift in DOS is 
observed when Mg is removed from MgB2. 

BeB2 could be expected to have a high owing to the lighter Be atoms which may 
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provide larger phonon frequencies while maintaining a similar electronic structure to that 
of MgB2. Even though BeB2 is isoelectronic with MgB2, the recent experimental study!! 
did not reveal any sign of superconductivity down to 5K. This negative finding makes it 
interesting to investigate the electronic structure of BeB2 in detail. The lattice constants for 
BeB2 obtained^ by averaging experimental data for the actual unit cell have been inferred 
to be a = 2.94 and c = 2.87 A. There are no experimental lattice parameters available for 
BeB2 with the AlB2-type structure and hence we have made structural optimization using 
the FPLAPW method. Table|I| shows that a is reduced only by about 1.6 %, whereas c is by 
about 8 % when Mg is replaced by Be. The anisotropic changes of the lattice parameters 
can be understood from the anisotropic bonding situation in MgB2. It looks as if the strong 
B {p — p) cr bonds within the planes prevent more appreciable changes in a whereas the 
weaker tt bonds along with the relatively weak ionic bonding between Be and B bring about 
a larger change in the c. Consistent with the earlier calculationS our calculated value of 
DOS at Ep for BeB2 is smaller than that of MgB2. The smaller N{Ep) value results from 
the broad nature of band structure of BeB2 compared with that of MgB2. Owing to the 
smaller volume our calculated bulk modulus for BeB2 becomes larger than that of MgB2 
(see Table 1). 

The total DOS curve for BeB2 (see Fig.^) shows almost free-electron-like metallic feature, 
with a N{Ef) value of 6.309 state Ry~^ f.u.~^ which in turn consistent with the paramag- 
netic behavior observed experimentally. The DOS curves of BeB2 and MgB2 are indeed 
very similar as expected due to their isoelectronic and postulated isostructural nature. As 
mentioned above the DOS features could lead one to a higher Tc for BeB2 than for MgB2 
(the molecular weight of BeB2 is lower than that of MgB2). Hirish predicted^ that the 
charge transfer from Be to B in BeB2 is less than that from Mg to B in MgB2 and that 
Ep is below the regime where superconductivity occurs. A closer inspection of the band 
structure (see Fig. § indicates that the key energy band, which we believe to be responsible 
for superconductivity, is broader in BeB2 than in MgB2 and it is also located well above Ep- 
Hence the calculation suggests that even if one stabilizes BeB2 in the AlB2-type structure 
one can not expect superconductivity. Our conclusion is consistent with the experimental 
observation in the sense that a recent studyi§ shows paramagnetic behavior down to 5 K. 

Kortus et a/.§ suggested that Ca doping should lead to an overall increase in DOS, and 
also provide an additional contribution to the electron-phonon coupling. We have therefore 
calculated DOS (using the optimized structural parameters) for CaB2 which shows (Fig.^ 
sharp features like those found in transition metal phases resulting from enhancement in 
volume compared with MgB2. Also, the calculated DOS at Ep is larger than that of MgB2 
indicating a possibility for superconductivity. The larger volume compared with BeB2 and 
MgB2 along with the weak B-B interaction make the bulk modulus for this material become 
smaller. The electronic structure of CaB2 (Fig.|^) shows that the key energy band is broader 
than that in MgB2. Also this doubly degenerate band is well below Ep at the F point 
and well above Ep at the A point. Hence, the calculations suggest that the probability for 
superconductivity in CaB2 with AlB2-type structure is low. 

The total DOS for SrB2 (see Fig.|D predicts a pseudogap feature at Ep (separating 
bonding from antibonding states with a negligible DOS at Ep). The electronic structure of 
SrB2 (Fig.^ shows nearly semimetallic feature. Earlier studiesS'll indicate that materials 
with Ep located at a pseudogap in DOS will have relatively high stability. (The situation 
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with all bonding orbitals filled and all antibonding orbitals empty implies extra contribution 
to stability.) Hence the calculations predict that SrB2 with AIB2 structure may be stabilized 
experimentally if the above criterion works. However, materials with Ep located in the 
pseudogap are not expected to become .uperconductingB and hence the present hnding 
suggests that SrB2 should be non-superconducting. Compared to BeB2 and MgB2, the top 
of the VB in CaB2 and SrB2 have large nonbonding B p states (see Fig.^). This will give 
negative contribution to the one-electron eigen-value sum for stability. This may be the 
reason why no stable CaB2 and SrB2 compounds are stabilised experimentally in the AIB2 
structure. 

The number of valence electrons per B, C atom is same for LiBC and MgB2. So, if 
electron per atom ratio is the decisive factor for the superconducting behavior of MgB2 one 
can expect superconductivity for LiBC. Hence, we have also performed electronic structure 
studies for LiBC. The calculated DOS (Fig.^) predicts insulating behavior with a band gap 
of L81eV. LiBC is an indirect bandgap insulator where the bandgap is between the top of 
VB in the F-K direction and the bottom of CB at the M point (Fig. The establishment 
of insulating behavior is consistent with the experimental observationS of a very small 
conductivity for LiBC. The fiat band feature present in the K-M and F-A directions of the 
BZ just below Ep around —0.23 eV (see Fig. ^) suggests that LiBC may be tuned to become 
superconducting upon hole doping. 

The calculated total DOS (Fig. H) for MgB2C2 shows a finite number of electrons at Ep 
(which has the same number of valence electrons per non-metal atom as MgB2) and predicts 
metallic behavior. This is further confirmed by the band structure (Fig. which shows that 
several bands cross the Fermi level. The most interesting aspect of Fig. ^ is that there is fiat 
band present in the vicinity of Ep (around the T— Y direction) similar to that in MgB2 and 
superconducting transition metal borocarbides (e,g, YNi2B2C, LuNi2B2C). Consequently we 
suggest superconductivity with relatively high Tc for MgB2C2. A more detailed analysis of 
the band structure shows that the narrow band in the T— Y direction near Ep is stemming 
from the C Pz electrons. The calculated density of states at Ep for MgB2C2 is 10.98 states 
Ry~^ f.u.~^ which is higher than for MgB2. Except for the structural dataS no other 
information on the physical properties of MgB2C2 are available experimentally which may 
be a promising candidate for high superconductivity. 

If it is the boron layers that are responsible for the superconductivity, one could expect 
superconductivity in MgB4. However, the calculated total DOS for MgB4 shows (Fig.|D 
only features which point towards insulating behavior and hence superconductivity is not 
expected for this material. 

B. Superconductivity 

Let us first look for similaritiies between MgB2 and other superconducting materials. 
There appears to be different correlations^ between superconducting transition temperature 
and the average electronegativity {■)]) for high Tc cuprates and conventional superconductors. 
Conventional superconductors take rj value between 1.3 and 1.9 whereas that for high Tc 
cuprates falls in the range 2.43 to 2.68.0 The i] value for MgB2 is 1.733 which places this 
material in the category of the conventional superconductors. However, using the proposed 
correlation^ between superconducting transition temperature and the value of electronic 
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specific heat coefficient (7) for various superconductors we find tliat 7 for MgB2 falls in the 
region of the high Tc cuprates (viz. MgB2 has a low value of N{Ep) as well as high Tc). In the 
conventional superconductors (Al5-type) with A3B composition, the mutually orthogonally 
permuted linear ..A-A-A.. chains are believed to be responsible for the superconductivity. 
In MgB2, the B atoms in the zig-zag chains are believed to be playing an important role 
for superconducting behavior. For higher Tc of superconducting intermetallics ButleiH has 
suggested a rule of togetherness which prescribes that (for a given crystal structure, electron- 
per-atom ratio and period) the electron-phonon coupling is enhanced when the transition 
metal atoms are brought closer together. It is interesting to note that the B atoms are 
brought close together by the strong covalent bonding in MgB2. Earlier studied^! show 
that /3-ZrNCl is a semiconductor with a band gap of ~3eV that upon Li intercalation 
(the electron donor element) becomes superconducting with Tc = 13 K.il This similarity is 
present for MgB2 also, that the electron donation from Mg to B2 leads to superconductivity. 
Ledbetter§ points out that the superconducting transition temperature in high Tc cuprates 
increases with increasing Debye temperature (Od). The high Tc in MgB2 may be associated 
with the large Od in this material and this criterion would place MgB2 with high Tc cuprates. 

In conventional superconductors, Tc increases with decreasing 6d, i.e. with lattice 
softening.i3 The calculated Od for MgB2 from the elastic constants is 1016 K which is found 
to be much higher value than in phonon mediated superconducting materials. However, the 
calculated value of 60 is comparable with the experimental values obtained from specific 
heat measurements (Table |I|). The presently derived is exceptionally large indicating 
that some novel mechanism is responsible for the high Tc in MgB2. Concerning the relation- 
ship of the purely phononic properties to superconductivity, owing to the large 9d, average 
phonon properties cannot reliably be used to estimate superconducting properties for this 
material. Using the concept that temperature-dependent electronic screening arising from 
narrow bands in the vicinity of Ep causes temperature-dependent phonon-mode frequencies 
(see Fig.^, one is lead to expect that softening of phonon modes does occur. Therefore, ex- 
perimental temperature dependent phonon spectra for MgB2 is required to establish whether 
softening of particular phonon modes is responsible for its large Tc. 

The BCS theory and its subsequent refinements based on Eliasberg equations show that 
high Tc in phonon-mediated superconductors is favored by high phonon frequency and a large 
DOS at Ep. High Tc materials commonly display interesting peculiarities in the phonon- 
dispersion curves, often in the form of dips (i.e. softening of phonons in well-defined regions 
of reciprocal space) which are particularly evident when phonon spectrum of a high-Tc 
material is compared with that of a similar low-Tc material (such as Nb vs. Mo, TaC vs. 
HfC). Assuming that this is the case for MgB2 the high Tc can be explained as follows. The 
superconducting transition temperature for the strongly coupled superconductors according 
to the McMillan's formula^! is, 



In Eqn. |1IIB| large Tc can be obtained when we have a large value for Op) (which MgB2 has) 
and the electron-phonon coupling constant A. The empirical value of the Coulomb coupling 
constant fi* for sp metals is 0.1. The McMillan- Hopfield expression for A = ^^m<u/2> ^ ' which 
enters in the exponent in the expression for Tc. Here, N{Ep) is smaller for MgB2 than for 
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conventional high Tc materials, < > is the averaged square of the electron-phonon matrix 
element, < c<j^ > is the averaged square of the phonon frequency and M is the mass of the 
ion involved. If the superconductivity occurs by phonon mediation, high in MgB2 can 
be explained as follows. Selected phonon modes may be strongly coupled to the electronic 
system and influence the magnitude of Tc to a greater extent than averagenhonon modes. In 
consisitent with the above view point the lattice dynamical calculation^iloEd'Ll reveal that 
the in-plane boron phonons near the zone-center are highly anharmanic with significant 
non-linear contribution to the electron-phonon coupling. If such a situation occurs, the 
softening of certain phonons in combination with the light mass of the constituents may 
lower the phonon contribution (M < cu^^ >) and in turn enhance the electron-phonon 
coupling constant and hence Tc. This may explain why MgB2 possesses high Tc despite the 
lower value of N{Ef). 

Now we will compare MgB2 with relevant superconducting material. Among the other 
transition metal diborides, superconductivity with Tc=9.5K has been observed in TaB2 
recently.Ell Our electronic structure studieB show that the Ep is located in a peak in the 
DOS profile in TaB2. This feature along with the large N{Ef) value of 12.92 states/ (Ry f.u.) 
will explain the large Tc in this material. As there is substantial B p states present at Ep in 
MgB2 similar to superconducting RNi2B2C, it may be worthwhile to compare these two cases. 
An experimental soft x-ray emission spectroscopy studylll on superconducting YNi2B2C and 
non-superconducting LaNi2B2C along with bandstructure calculationEl indicate that the 
superconductivity appears only when the broad B p bands are located at Ep- The main 
difference between the electrons involved in transport properties of RNi2B2C and MgB2 is 
that the former has a remarkable DOS peak at Ep dominated by Ni 3d states with almost 
equal proportions of all five Ni 3d states and also involves some rare-earth d and B, C 
sp admixture whereas MgB2 does not show any peak feature at Ep. From the accurate 
analysis of measure specific heat over a wide temperature range the estimated 7 for 
MgB2 is estimated as 5.5mJ/mol K^. From our calculated N{Ep) value we have derived 
the electronic contribution to the specific-heat coefficient without electron-phonon mass 
enhancement yield a value of 1.73mJ/mol K^. From this value along with the experimental 
lexp we have estimated the value of electron-phonon coupling constant using the relation 
lexp = 7t/i(l + A) which gave A = 2.17. The large A value indicates MgB2 is a strongly 
coupled superconductor. This strong coupling as well as the low mass of B can explain the 
high Tc in MgB2. 

Using the calculated plasma frequency (flpx) in the ab plane (7.13 eV) along with N{Ep) 
= 0.71 states/(eV f.u.) the in-plane Fermi velocity up^ has been obtained as 



fi2 

- 8.9 X 10^ cm/s. 



\ 47re'^N{Ep 



Using the superconducting gap A=3.53A;bTc/2 = 6meV we have calculated the coherence 
length (^) and the field penetration depth (A) 

^ = = 307 A; A = = 276A. 

ttA ilpx 

The experimental upper critical field, Hc2(T), thermodynamic critical field, Hc{T), and 
critical current, Jc, indicate that MgB2 is a type-II superconductor.0 
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Now we will try to understand the effect of pressure on the superconductivity in MgB2. 
Addition of Mg to the boron sublattice gives 15% shortening of the c axis compared to that 
in graphite.0 This chemical pressure brings the B semimetal to superconducting state. As 
the transport properties of MgB2 is from the holes, it is suggested that the superconductivitv 
may be understood within the formalism developed for high-Tc cuprate superconductors.!^ 
This predicts a positive pressure coefficient for Tc as a result of the decreasing intraplane B-B 
distance with increasing pressure. However, the Tc is experimentally0 found to decrease with 
pressure at a rate of — 1.11 K/GPa. The RVB theory,0 on the other hand, predicts a decrease 
in Tc with increasing chemical pressure, consistent with the experimental findingsHil00 
Linear response calculations^ show that the B-bond stretching modes have unusually strong 
coupling to electrons close to Ep at the top of the bonding quasi-two-dimensional B cr bonds. 
Therefore, when the compressibility is larger within the ab plane than along c, one can expect 
a large variation in the superconductivity with pressure (because the flat bands with p^^y 
characters broaden faster by compression). Our calculated elastic property for MgB2 shows 
large anisotropy with easy compression along c. Hence, the key energy band (which is 
sensitive to the B bond stretching) will broaden slowly with pressure and Tc will slowly 
decrease with increasing pressure, consistent with experimental results. 

From the FPLAPW calculations we have estimated the electric field gradient (V^^) at the 
Mg and B site as —0.249 xlO^^ and 2.047 x lO^^V/m^, respectively. Using the calculated 
Vzz along with the nuclear quadrapole moment for ^^B (0.037 we can calculate the NMR 
quadrapole coupling frequency (z/g) by means of the relation 



2hl{21 -1) 

using I = 3/2 as the nuclear spin quantum number for ^^B. This gave Ug = 915 kHz that is 
in good agreement with 828 kHZ obtained by Gerashenko et al^ and 835 ± 5 kHz obtained 
by Jung et al. from first order quadrapole perturbed NMR spectrum.lll 



C. Chemical bonding in MgB2 



Similar to the C-C distances in the graphite structure, the distance between the boron 
planes in MgB2 is about twice the intraplanar B-B distance and hence the B-B bonding 
is strongly anisotropic. A more quantitative assessment of the bonding situation in MgB2 
can be obtained from the partial DOS (Fig.^) which demonstrates that the B s states 
are hybridized with the B p state in VB. This shows strongly bonded sp^ hybrids in the 
ah plane. The Mg s electrons contribute very little to the VB and are mainly reflected 
in the unoccupied state. Hence Mg donates electrons to the boron layers. VB carries 
predominantly B 2p character formed from two distinct sets [a (jix,y) and tt {pz)] of bonds. 
The B 2s electrons are well localized and their contribution at Ep is minor. From the 
orbital projected DOS (Fig.|^) it is clear that B p^ and Py characters are mainly dominating 
at Ep. This suggests that the p — p a bonding between the boron atoms has a significant 
influence on superconductivity. It is worth to recall that the recently found@ superconductor 
Nao.29HfNCl (Tc ~ 25 K) has N p^ and Py characters at E'^^.S B pz states are present in a 
wide energy range and dominating at the bottom of CB. 
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The simplest way to investigate the bonding situation between two interacting atoms in 
sohd is to inspect the complete COHP between them, taking all valence orbitals into account. 
In the upper panel of Fig|^ shows COHP for the B-B and Mg-B bonds. An interesting aspect 
of this illustration is that VB is filled up with bonding orbitals (negative value of COHP) 
and the antibonding orbitals are some ~3eV above Ep. Bonding state electrons from both 
B-B and Mg-B bonds are found at Ep- The B s — s bonding states are found mainly at 
the bottom of VB around — 8.5eV. The B p — p a bonding states dominate at the top of 
the VB region around — 2eV. Note that the COHP values for the Mg-B bonds are much 
smaller than for the B-B bonds indicating that the B-B bond is much stronger than the 
Mg-B bonds. In order to quantify the bonding interactions in MgB2 we have integrated 
the COHP curve and obtained the values —4.36, —1.23 and — 0.307eV for B-B, Mg-B and 
Mg-Mg bonds respectively. This further confirms that the B-B bonds are the strongest in 
MgB2 which in turn consistent with the derived elastic properties (see Sec. [HI D| ). 

In order to further illustrate the bonding situation in MgB2 a charge density plot for the 
(110) plane is shown in Fig.|^. This illustration shows a low electron accumulation between 
Mg and B as well as a very low electron population at the Mg site (much lower than for 
a neutral Mg atom). These findings are a clear indication of ionic bonding between Mg 
and B. The large electron accumulation between the B atoms and their strongly aspherical 
character indicate strong covalent interaction between the B atoms as also found by our 
examination of partial DOS and COHP. The more or less homogeneous charge distribution 
between the Mg atoms suggest an appreciable degree of metallic bonding between them. i.e. 
apart from strong ionic and covalent bonding in MgB2 the band structure shows features 
similar to the sp metals. Hence, MgB2 is a typical example of a mixed bonded solid. 

D. Elastic properties 

Structural parameters for BeB2, CaB2, SrB2 are not available experimentally and we 
have therefore made structural optimization for these diborides along with the MgB2 by total 
energy minimization. The calculated total energy variation as a function of c/a and volume 
for MgB2 is shown in Fig. ^ and §3 respectively. The optimized structural parameters along 
with the bulk modulus, its pressure derivative and DOS at Ep for all these diborides are 
given in Table |. The calculated equilibrium volume for MgB2 is found to be in excellent 
agreement with the experimental value,@ and the corresponding c/a value is only 0.42% 
smaller than the experimental value (see Fig. ^ and Table|). These results indicate that 
density functional theory works well for this material. The calculated N{Ep value for MgB2 
is larger than that for other diborides considered in the present study and this may be one 
of the reasons for the superconductivity in MgB2. 

By fitting the total-energy vs. volume curve in Fig.^ to the universal equation of state 
we have obtained the bulk modulus (Bq) and its pressure derivative (B'q) as 1.50 Mbar and 
3.47, respectively. The derived Bq is found to be in excellent agreement with the value 
1.51 Mbar obtained from high resolution x-ray powder diffraction measurement J^l and also 
in good agreement with other results listed in Table|l| Furthermore Bo = 1.5 Mbar obtained 
in this way is in good agreement with the value 1.509 Mbar estimated from the calculated 
single-crystal elastic constants by means of the relation 
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^ _ 2(cii + C12) + 4ci3 + C33 
9 

The electron per atom ratio of TiB2 is same as that of MgB2 and hence this material should 
also be considered. In TiB2 Ep is located in a pseudogap and hence N{Ef is small resulting 
in non-observation of superconductivity. The calculated bulk modulus for TiB2 is larger 
than that for diborides considered in the present study. This is an effect of the filling of the 
bonding states in the bands and this aspect has been discussed in Ref. 

From various distortions we have calculated all the five single- crystal elastic constants 
for MgB2 (see Table |ITTD. Unfortunately there are no experimental elastic constants available 
since a suitable single crystal of MgB2 has so far not been obtained. We have therefore made 
comparison with the data for TiB2 available experimentallyEil. Despite the layered crystal 
structure of MgB2 the high pressure total energy studies of Loa and SyassenEHl suggested 
isotropic compressibility and they concluded that the intra- and the interlayer bonding are 
of similar strength. Another high pressure study (upto 8 Gpa) by Vogt et alx3 concluded that 
there are small anisotropics in the mechanical properties. The isothermal compressibility 
measurements of MgB2 by synchrotron x-ray diffraction revealed^ a stiff tightly-packed 
incompressible nature with only moderate anisotropy between intra- and interlayer bonds. 
The neutron diffraction measurements by Jorgensen et a/.il at high pressures concluded 
with highly anisotropic mechanical properties for MgB2. As the experimental results are 
mutually inconsistent, theoretical studies of the compressibility may be helpful in resolving 
the ambiguities. 

From the calculated single- crystal elastic constants we have derived the bulk moduli along 
the crystallographic directions using the relations given in Sec. || and the results are listed in 
Table pn]. The compressibility studyS by synchrotron radiation show isothermal interlayer 
compressibility, dlnc/dP at zero pressure is 1.4 times the inplane compressibility, dlna/dP. 
However the high pressure neutron diffraction studieJ^l show that dlnc/dP=1.64 dlna/dP 
which is much closer to our calculated relationship dlnc/dP=1.79 dlna/dP (see Table pT|. 
The most recent high pressure measurement^ by hydrostatic pressure upto 15 GPa shows a 
large anisotropy with the relationship dlnc/dP=1.875 dlna/dP. The higher compressibility 
along c than along a can be understood as follows. There is strong B p^ y — p^^y covalent 
hybridization along a in MgB2 and hence the bulk modulus along a {b) is large. There is 
significant ionic contribution to the bonding between Mg and B along c. Usually an ionic 
bond is weaker than a covalent bond and hence the bulk modulus is smaller along c than 
along a. The large anisotropy in the compressibility is also consistent with the fact@ that 
the thermal expansion of c is about twice as that of a. Moreover, substitutionlli of Al for 
Mg site decreases c at a rate approximately twice that of a indicating the anisotropic nature 
of bonding. 

It is possible to visualize the anisotropy in the elastic properties from the curvature of 
the total energy with respect to length changes in an arbitrary direction. From the elastic 
compliance constants (sjj) it is possible to derive the directional bulk modulus K, using the 
following relation:ii 

^ = {Su + S12 + S13) - (Sii + S12 - Si3 - S33)/3 (4) 
A 

where I3 is the direction cosine. Thus obtained directional dependent bulk modulus is shown 
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in Fig. |a. A useful surface construction is one that shows the directional dependence of 
Young's modulus (E), which for hexagonal symmetry can be defined as 

i = (1 - ll)hn + 4^33 + llil - ll)i2s,s + su) (5) 

Thus derived directional-dependent Young's modulus (Fig. |b) also shows large 
anisotropy, E along a being about 65% larger than along c. The anisotropic nature of 
bonding behavior reflected in the elastic properties is consistent with the charge density 
analysis. The marked anisotropic compressibility of MgB2 will lead to different pressure 
effects on different phonon modes and is also more likely to lead to pressure induced changes 
in the electronic structure at Ep. This information is valuable in testing the predictions of 
competing mo delffl for the mechanism of superconductivity. 

The anisotropy in the plastic properties of materials can be studied from the directional 
dependent shear stress and the amount of shear.il So, it is interesting to study the directional 
dependence of shear in MgB2. We have calculated the directional dependent shear modulus 
(G) from the elastic constants using the relation 

^ = S44 + (Sll - Sl2 - ^)(1 - ^3) + 2(S11 + S33 - 2Si3 - S44)(l - ^3)^3 (6) 

The directional dependent shear modulus (see Fig.|^) also shows large anisotropy and the 
G along the a axis is around 42% higher than that along c. It should be noted that large 
shear value is present in between the basal plane and perpendicular to basal plane. From 
the calculated directional dependent bulk modulus and the shear modulus one can calculate 
the directional dependence of sound velocity and hence the characteristic temperature of the 
material. The calculated characteristic temperature given in Fig. ^ shows anisotropic nature 
due to the anisotropy in the bonding behavior. The calculated average shear modulus (Gav), 
elastic-wave velocities (z/) and 6d obtained from our single crystal elastic constants are listed 
in Table m The calculated On for MgB2 is much higher than the experimental values and the 
higher value of do indicates that higher phonon frequencies play an important role for the 
high Tc of this material. The existence of phonon modes at very high energies as 97meV in 
MgB2 is experimentally identified by neutron-inelastic-scattering measurements.ElS One of 
the reasons for the discrepancy between the 6d obtained from the elastic constants and the 
specific heat measurements is as follows. For example, the experimental electrical resistivity 
and specific heat measurements on TiSi2 gave 9d = 560-664 K, whereas elastic constants 
measurements gave 6d = TOOK (that is in close agreement with Op, = 722 K obtained from 
our theoretical elastic constant calculationsH). So, one can expect good agreement when 
comparing our On with that obtained from elastic constants measurement. 



E. Optical properties 

Optical properties studies are of fundamental importance, since these involve not only the 
occupied and unoccupied parts of the electronic structure but also carries information on the 
character of the bands. In order to elucidate the anisotropy in the optical properties of MgB2 



the calculated imaginary parts of dielectric tensor for E\\a and E\\c are shown in Fig. 10 



The important feature conveyed by the interband transition shown in Fig. 10 is that there is 
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negligible contribution to €2(00) in E\\c below 3.8 eV. The interesting aspect for MgB2 is that 
even though the interband contribution to the optical spectra is highly anisotropic the calcu- 
lated intraband contribution show nearly isotropic behavior. This originates from the close 
values of the calculated plasma frequencies; 7.13 and 6.72 eV for in-plane and perpendicular 
to the plane, respectively. These values are in excellent agreement with Up^x= 7.02 eV and 
^p,z= 6.68 eV obtained by full-potential LMTO calculation.0. The intraband contribution 
to the optical dielectric tensor have been calculated similar to our earlier study.ii The calcu- 
lated e2(w) spectra which include both inter- and intraband contributions are shown in the 
upper panel of fig.|TU[ From these €2(0;) spectra one can derive all linear optical properties. 
Unfortunately, there are no experimental spectra available for MgB2. 

IV. SUMMARY 

The present study shows a common origin between the superconductivity in rare earth 
transition metal borocarbides and MgB2 that there is a flat band present in the vicinity 
of Fermi level and also the B atoms are primarily involved in electron-phonon coupling in 
these materials. Owing to the lack of single crystals the anisotropy in physical properties 
of MgB2 are not studied experimentally. We predicted large anisotropics in the optical and 
mechanical properties. From detailed electronic structure studies we have arrived at the 
following conclusions. 

1. The bonding behavior in MgB2 has been explained by analyses of site, angular momentum 
and orbital projected density of states as well as charge density and crystal overlap Hamil- 
tonian population. These analyses establish a mixed bonding behavior with ionic bonding 
between Mg and B, covalent bonding between B atoms and metallic bonding between Mg 
and Mg like that in sp metals. 

2. We identified a large anisotropy in the mechanical properties of MgB2 from our calcu- 
lated elastic constants, consistent with the anisotropy in the bonding behavior and high 
pressure neutron diffraction measurements. Consequently, pressure can influence the bands 
in different directions of the BZ in an unusually different manner and hence the physical 
properties. 

3. Two degenerate flat bands have been identified near Ep in the F-A direction of BZ. These 
degenerate B p^^y bands are considered as the key in the realization of the high temperature 
superconductivity of MgB2. 

4. The role of Mg in MgB2 is to donate electrons to B atoms and hence to shift Ep such 
that it lies very closer to the flat band, the feature which we believe important for super- 
conductivity. 

5. We found similarity in the electronic structures of MgB2 LiBC and MgB2C2. Therefore 
superconductivity is expected in MgB2C2 and hole doped LiBC. On the other hand the 
electronic structure of BeB2, CaB2, SrB2 and MgB4 suggest a low probability for supercon- 
ductivity. 
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6. Because of the nearly flat bands in the vicinity of Ep (which will introduce temperature 
depending electronic screening of the phonon-modc frequencies), we expect temperature- 
dependent phonon-mode softening or large anisotropy in the phonon modes. 

7. The calculated boron NMR frequency is found to be in very good agreement with the 
experimental studies. 

8. As our calculated Debye temperature for MgB2 is much larger than for other supercon- 
ducting intermetaUics, we believe that selected phonon modes can be strongly coupled to 
the electronic system and thus influence the magnitude of Tc to a greater extent than the 
average phonon correlations may indicate. 

9. If the rigid band- filling approximation works and also the electrons are responsible for 
superconductivity in MgB2, our calculation suggests that a doping of 0.32 electrons will 
enhance Tc. However, the experimental observation of reduction in by electron doping in 
MgB2 indicate that the holes are responsible for superconductivity. 
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FIGURES 



FIG. 1. Crystal structures of MgB2, MgB2C2, LiBC and MgB4. Legends to the different kinds 
of atoms are given on the illustration. 

FIG. 2. Band structure of MgB2 and hypothetical B2 in the AlB2-type framework. For the 
B2 substructure we have used the lattice parameters of MgB2 . 

FIG. 3. Lower panel: Total and site projected DOS for MgB2 and hypothetical B2 obtained 
from the FPLMTO method. Middle panel: Angular momentum and orbital projected DOS for 
MgB2 obtained from the FPLAPW method Upper panel: Crystal overlap Hamiltonian population 
(COHP) between Mg-Mg, Mg-B and B-B obtained from the TBLMTO method. 

FIG. 4. Calculated total DOS for BeB2, CaB2, SrB2 obtained from the FPLAPW method 
and that for MgB4, LiBC and MgB2C2 obtained from the TBLMTO method. For BeB2, CaB2, 
and SrB2 we have assumed the AlB2-type structure and the optimized structural parameters. 

FIG. 5. The band structure of BeB2, CaB2 and SrB2 with AlB2-type structure as obtained 
using the optimized lattice parameters. 

FIG. 6. The band structure of LiBC and MgB2C2 as obtained with the TBLMTO method 
using experimental lattice parameters. 

FIG. 7. Valence charge density in the (110) plane for MgB2. Corner atoms are Mg and the 
others are B. Both plots have 35 contours between and 0.25 electrons/a.u.^. 

FIG. 8. The variation of the total energy with unit-cell volume (lower panel) and c/a (upper 
panel) for MgB2 (AlB2-type structure) as obtained from the FPLMTO calculations. 

FIG. 9. Calculated directional-dependent (a) bulk, (b) Young's, (c) shear moduli and the 
characteristic temperature (d) for MgB2 as obtained from the calculated single crystal elastic 
constants. 

FIG. 10. The calculated optical dielectric tensor for MgB2 obtained from with and without 
intraband contribution obtained from the FPLMTO calculation. 
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TABLES 



TABLE L Calculated lattice parameters (a and c are in A), c/a ratio, bulk modulus (Bq in 
Mbar), its pressure derivative (Bg) and density of state at the Fermi level [N{Ep) in states Ry~^ 
f.u.~^] for AlB2-type compounds 



Compound 


a 


c 


c/a 


Bo 


B'o 


N{Ef) 


BeB2 


2.886 


3.088 


1.027 


1.93 


3.47 


6.309 


MgBa 


3.080 


3.532 


1.147 


1.50 
1.20±0.05^ 

1.39^ 

1.5r 

1.47"^ 
1.40±0.06*^ 


3.50 


9.98 


CaB2 


3.397 


4.019 


1.183 


1.34 


3.42 


10.837 


SrBs 


3.456 


4.193 


1.213 


1.05 


2.28 


4.732 


TiBa 


3.070 


3.262 


1.060 


2.13 


2.10 


4.27 



'^Prom synchrotron XRD by Prassides et a/.0. 

^ High pressure XRD measurement by Vogt et aZ.il. 

Calculated from FPLAPW method by Vogt et a/.ii. 

Calculated from pseudopotential method by Bohnen et alx 
^ Calculated from FPLAPW method by Lao and Syassenll. 



TABLE II. Average shear modulus {Gav in Mbar), longitudinal and transverse elastic wave 
velocity {vi , i/^ in m/s) and Debye temperature {Qd in K) obtained from single crystal elastic 
constant (Od in K) 

Compound Gav i^i i^m ©^(calc.) 0z)(exp.) 

MgB2 1.146 10612 7276 1016 750±30^ 

746^ 
800^= 
920^ 

TiB2(exp.)<' 2.407 7061 5105 743 

^ From specific heat measurements by Bud'ko et a/.0. 
^ From specific heat measurements by Walti et al^. 
^ From specific heat measurements by Kremer et a/.i3. 

From specific heat measurements by Wang et aZ.il. 
^Calculated from the single crystal elastic constants of Spoor et aZ-S 
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TABLE III. The single crystal elastic constants {cij in Mbar) and bulk modulus values along 
a and c {Ba and Be in Mbar) 



Compound 


Cll 


Cl2 


Cl3 


C33 


C44 


Ba 


Be 


MgBa 


4.380 


0.430 


0.329 


2.640 


0.802 


5.406 


3.006 














4.1±0.2'^ 


2.92±0.12 














6.25'' 


3.33 


TiB2(exp.)" 


6.6 


0.48 


0.93 


4.32 


2.6 


8.512 


5.527 



'^From synchrotron XRD by Prassides et al. 

''Prom hydrostatic high pressure synchrotron XRD by Goncharov et all 
'^Calculated from the single crystal elastic constants of Spoor et a/.@. 
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r KM rA LH A 

LiBC 



r z T Y r 

MgB.C, 



